Esx1 is an X-chromosome-imprinted regulator of placental formation and fetal growth

Placental formation and genomic imprinting are two important features of embryonic development in placental mammals. Genetic studies have demonstrated that imprinted genes play a prominent role in regulating placental formation. In marsupials, mice and humans, the paternally derived X chromosome is preferentially inactivated in the placental tissues of female embryos. This special form of genomic imprinting may have evolved under the same selective forces as autosomal imprinted genes. This chromosomal imprinting phenomenon predicts the existence of maternally expressed X-linked genes that regulate placental development.^ In this study, an X-linked homeobox gene, designated Esx1 has been isolated. During embryogenesis, Esx1 was expressed in a subset of placental tissues and regulates formation of the chorioallantoic placenta. Esx1 acted as an imprinted gene. Heterozygous female mice that inherit an Esx1-null allele from their father developed normally. However, heterozygous females that inherit the Esx1 mutation from their mother were born 20% smaller than normal and had an identical phenotype to hemizygous mutant males and homozygous mutant females. Surprisingly, although Esx1 mutant embryos were initially comparable in size to wild-type controls at 13.5 days post coitum (E13.5) their placentas were significantly larger (51% heavier than controls). Defects in the morphogenesis of the labyrinthine layer were observed as early as E11.5. Subsequently, vascularization abnormalities developed at the maternal-fetal interface, causing fetal growth retardation. These results identify Esx1 as the first essential X-chromosome-imprinted regulator of placental development that influences fetal growth and may have important implications in understanding human placental insufficiency syndromes such as intrauterine growth retardation (IUGR). ^

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Genomic screening identifies novel linkages and provides further evidence for a role of MYH9 in nonsyndromic cleft lip and palate.

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth anomaly that requires prolonged multidisciplinary rehabilitation. Although variation in several genes has been identified as contributing to NSCLP, most of the genetic susceptibility loci have yet to be defined. To identify additional contributory genes, a high-throughput genomic scan was performed using the Illumina Linkage IVb Panel platform. We genotyped 6008 SNPs in nine non-Hispanic white NSCLP multiplex families and a single large African-American NSCLP multiplex family. Fourteen chromosomal regions were identified with LOD>1.5, including six regions not previously reported. Analysis of the data from the African-American and non-Hispanic white families revealed two likely chromosomal regions: 8q21.3-24.12 and 22q12.2-12.3 with LOD scores of 2.98 and 2.66, respectively. On the basis of biological function, syndecan 2 (SDC2) and growth differentiation factor 6 (GDF6) in 8q21.3-24.12 and myosin heavy-chain 9, non-muscle (MYH9) in 22q12.2-12.3 were selected as candidate genes. Association analyses from these genes yielded marginally significant P-values for SNPs in SDC2 and GDF6 (0.01

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

A conserved motif at the 3$\sp\prime$ end of genomic RNA of mouse hepatitis virus required for host protein binding and viral RNA replication

The initial step in coronavirus-mouse hepatitis virus (MHV) replication is the synthesis of negative strand RNA from a positive strand genomic RNA template. Our approach to studying MHV RNA replication is to identify the cis-acting signals for RNA synthesis and the protein(s) which recognizes these signals at the 3$\sp\prime$ end of genomic RNA of MHV. To determine whether host cellular and/or virus-specific proteins interact with the 3$\sp\prime$ end of the coronavirus genome, an RNase T$\sb1$ protection/gel mobility shift electrophoresis assay was used to examine cytoplasmic extracts from either mock- or MHV-JHM-infected 17Cl-1 murine cells for the ability to form complexes with defined regions of the genomic RNA. A conserved 11 nucleotide sequence UGAAUGAAGUU at nucleotide positions 36 to 26 from the 3$\sp\prime$ end of genomic RNA was identified to be responsible for the specific binding of host proteins, by using a series of RNA probes with deletions and mutations in this region. The RNA probe containing the 11 nucleotide sequence bound approximately four host cellular proteins with a highly labeled 120 kDa and three minor species with sizes of 103, 81 and 55 kDa, assayed by UV-induced covalent cross-linking. Mutation of the 11 nucleotide motif strongly inhibited cellular protein binding, and decreased the amount of the 103 and 81 kDa proteins in the complex to undetectable levels and strongly reduced the binding of the 120 kDa protein. Less extensive mutations within this 11 nucleotide motif resulted in variable decreases in RNA-protein complex formation depending on each probe tested. The RNA-protein complexes observed with cytoplasmic extracts from MHV-JHM-infected cells in both RNase protection/gel mobility shift and UV cross-linking assays were indistinguishable to those observed with extracts from uninfected cells.^ To investigate the possible role of this 3$\sp\prime$ protein binding element in viral RNA replication in vivo, defective interfering RNA molecules with complete or partial mutations of the 11 nucleotide conserved sequence were transcribed in vitro, transfected to host 17Cl-1 cells in the presence of helper virus MHV-JHM and analyzed by agarose gel electrophoresis, competitive RT-PCR and direct sequencing of the RT-PCR products. Both negative strand synthesis and positive strand replication of DI RNA were affected by mutation that disrupts RNA-protein complex formation, even though the 11 mutated nucleotides were converted to wild type sequence, presumably by recombination with helper virus. Kinetic analysis indicated that recombination between DI RNA and helper virus occurred 5.5 to 7.5 hours post infection when replication of positive strand DI RNA was barely observed. Replication of positive strand DI RNAs carrying partial mutations within the 11 nucleotide motif was dependent upon recombination events after transfection. Replication was strongly inhibited when reversion to wild type sequence did not occur, and after recombination, reached similar levels as wild type DI RNA. A DI RNA with mutation upstream of the protein binding motif replicated as efficiently as wild type without undergoing recombination. Thus the conserved 11 nucleotide host protein binding motif appears to play an important role in viral RNA replication. ^

TUMOR PROGRESSION, METASTATIC POTENTIAL AND GENOMIC INSTABILITY

To investigate the hypothesis that increased malignant potential correlates with increased levels of genetic instability, the following parameters of instability were measured: (1) spontaneous mutation rates for ouabain resistance in murine cell lines of different malignant potentials, (2) the background prevalence of 6-thioguanine (6-TG) resistance in clone 4 (highly metastatic) and clone 19 (poorly metastatic) of the K1735 murine melanoma, (3) the prevalence of ouabain resistant variants in three murine cell lines and their variants after exposure to the mutagen MNNG, (4) the rate of generation of major karyotypic abnormalities in B16 F1 (poorly metastatic) and B16 F10 (highly metastatic) murine melanoma, and (5) analysis of the G-banded karyotypes of cloned B16 F1 and B16 F10 melanoma.^ No correlation of increased spontaneous mutation rates with increased malignant potential was found in repeated experiments with three murine cell lines and their variants of different malignant potential. The background prevalence of g-TG resistance was not significantly different for the poorly and highly metastatic clones of K1735 melanoma. The studies with MNNG-induced mutation showed no increased sensitivity of the highly metastatic variants of the three murine cell lines to mutagenesis. Neither did the rate of generation of major karyotypic abnormalities correlate with malignant potential. However, certain karyotypic differences were demonstrated after G-banding of the B16 F1 and F10 melanomas.^ One hypothesis which is consistent with these results is that the rate of generation of genetic abnormalities need not be strongly related to the degree of malignant potential. An increased prevalence of genetic changes may merely reflect the accumulation of abnormalities while their rate of production remains constant. The presence of specific nonrandom changes likely is the main determinant of malignant potential rather than the rate of production of random changes. ^

A FUNCTIONAL VARIANT OF HLA-DR1 DELINEATED BY T-LYMPHOCYTE CLONE REACTIVITY, PROTEIN CONFORMATION, AND GENOMIC RESTRICTION SITES

This research characterized a serologically indistinguishable form of HLA-DR1 that: (1) cannot stimulate some DR1-restricted or specific T-lymphocyte clones; (2) displays an unusual electrophoretic pattern on two dimensional gels; and (3) is marked by a polymorphic restriction site of the alpha gene. Inefficient stimulation of some DR1-restricted clones was a property of DR1$\sp{+}$ cells that shared HLA-B14 on the same haplotype and/or were carriers of 21-hydroxylase (21-OH) deficiency. Nonclassical 21-OH deficiency frequently demonstrates genetic linkage with HLA-B14;DR1 haplotypes and associates with duplications of C4B and one 21-OH gene. Cells having both stimulatory (DR1$\sb{\rm n}$) and nonstimulatory (DR1$\sb{\rm x}$) parental haplotypes did not mediate proliferation of these clones. However, heterozygous DR1$\sb{\rm x}$, 2 and DR1$\sb{\rm x}$, 7 cells were efficient stimulators of DR2 and DR7 specific clones, respectively, suggesting that a trans acting factor may modify DR1 alleles or products to yield a dominant DR1$\sb{\rm x}$ phenotype. Incompetent stimulator populations did not secrete an intercellular soluble or contact dependent suppressor factor nor did they express interleukin-2 receptors competing for T-cell growth factors. Two dimensional gel analysis of anti-DR immunoprecipitates revealed, in addition to normal DR$\alpha$ and DR$\beta$ chains, a 50kD species from DR1$\sb{\rm x}$ but not from the majority of DR1$\sb{\rm n}$ or non-DR1 cells. The 50kD structure was stable under reducing conditions in SDS and urea, had antigenic homology with DR, and dissociated after boiling into 34kD and 28kD peptide chains apparently identical with DR$\alpha$ and DR$\beta$ as shown by limited digest peptide maps. N-linked glycosylation and sialation of DRgp50 appeared to be unchanged from normal DR$\alpha$ and DR$\beta$. Bg1II digestion and $DR\alpha$ probing of DR1$\sb{\rm x}$ genomic DNA revealed a 4.5kb fragment while DR1$\sb{\rm n}$ DNA yielded 3.8 and 0.76kb fragments; all restriction sites mapped to the 3$\sp\prime$ untranslated region of $DR\alpha$. Collectively, these data suggest that DRgp50 represents a novel combinatorial association between constitutive chains of DR that may interfere with or compete for normal T cell receptor recognition of DR1 as both an alloantigen and restricting element. Furthermore, extensive chromosomal abnormalities previously mapped to the class III region of B14;DR1 haplotypes may extend into the adjacent class II region with consequent intrusion on immune function. ^

THE TRANSLATION PRODUCTS OF MOLONEY MURINE SARCOMA VIRUS 124 GENOMIC RNA

Murine sarcoma viruses constitute a class of replication-defective retroviruses. Cellular transformation may be induced by these viruses in vitro; whereas, fibrosarcomas may result in animals infected with them in vivo (Tooze, 1973; Bishop, 1978). Hybridization studies suggest that murine sarcoma viruses arose by recombination between nondefective murine leukemia virus sequences and certain cellular sequences present in uninfected mouse cells (Hu et al., 1977). A specific gene product, however, has not been implicated in murine sarcoma virus transformation.^ One line of murine sarcoma virus-producing cells, Mo-MuSV-clone 124, (Ball et al., 1973), was studied biochemically because it mainly produces the sarcoma virus as a pseudotype packaged with helper murine leukemia virus proteins. The sarcoma viral RNA was translated in a sophisticated cell-free protein synthesizing system (Murphy and Arlinghaus, 1978). The translation products were analyzed by a number of techniques, including electrophoresis in denaturing gels of SDS polyacrylamide, immunoprecipitation, and peptide mapping. The major products of the total RNA purified from the virus preparation were shown to have molecular weights of about 63,000 (P63('gag)), 42,000 (P42), 40,000 (P40), 38,000 (P38), and 23,000 (P23). The size class of mRNA coding for each of the cell-free products was estimated using a poly(A) selection technique and sucrose gradient fractionation. These analyses were used to localize the coding information related to each of the in vitro synthesized cell-free products within the sarcoma virus genome.^ The major findings of these studies were: (1) the 5' half of the sarcoma viral RNA codes for the 63,000 dalton polypeptide and 42,000 - 38,000 dalton polypeptides derived from the "gag" gene; and (2) the 3' half of the sarcoma viral RNA codes for a 38,000 dalton polypeptide and possibly derived from the cellular acquired sequences. ^

THE INHIBITORY EFFECTS OF VSV INFECTION ON MULV PRODUCTION OCCUR BETWEEN SYNTHESIS AND RELEASE OF MULV GENOMIC RNA

The inhibitory effects of VSV infection on MuLV production were investigated using the VSV temperature-sensitive mutants t1B17(I & V), tsT1026(I), tsG22(II), and ts052(II). At the permissive temperature, all four mutants suppressed the release of virion-associated MuLV gRNA by approximately 98% within 0.5 to 2.5 hr post infection. At the restrictive temperature and in the absence of cell killing, infection with t1B17(I & V) inhibited the release of MuLV gRNA, while tsT1026(I) and tsG22(II) did not. In contrast, ts052(II) inhibited the release of MuLV gRNA and induced cell killing. During the same time period and at either temperature, all four mutants did not suppress either MuLV-associated protein release or intracellular MuLV sRNA synthesis. These results indicate that VSV inhibits MULV gRNA release at a level somewhere between the synthesis and release of newly synthesized gRNA.^

Non B-DNA structures and genomic instability associated with myotonic dystrophy type 2 (CCTG).(CAGG) repeats

There are many diseases associated with the expansion of DNA repeats in humans. Myotonic dystrophy type 2 is one of such diseases, characterized by expansions of a (CCTG)•(CAGG) repeat tract in intron 1 of zinc finger protein 9 (ZNF9) in chromosome 3q21.3. The DM2 repeat tract contains a flanking region 5' to the tract that consists of a polymorphic repetitive sequence (TG)14-25(TCTG)4-11(CCTG) n. The (CCTG)•(CAGG) repeat is typically 11-26 repeats in persons without the disease, but can expand up to 11,000 repeats in affected individuals, which is the largest expansion seen in DNA repeat diseases to date. This DNA tract remains one of the least characterized disease-associated DNA repeats, and mechanisms causing the repeat expansion in humans have yet to be elucidated. Alternative, non B-DNA structures formed by the expanded repeats are typical in DNA repeat expansion diseases. These sequences may promote instability of the repeat tracts. I determined that slipped strand structure formation occurs for (CCTG)•(CAGG) repeats at a length of 42 or more. In addition, Z-DNA structure forms in the flanking human sequence adjacent to the (CCTG)•(CAGG) repeat tract. I have also performed genetic assays in E. coli cells and results indicate that the (CCTG)•(CAGG) repeats are more similar to the highly unstable (CTG)•(CAG) repeat tracts seen in Huntington's disease and myotonic dystrophy type 1, than to those of the more stable (ATTCT)•(AGAAT) repeat tracts of spinocerebellar ataxia type 10. This instability, however, is RecA-independent in the (CCTG)•(CAGG) and (ATTCT)•(AGAAT) repeats, whereas the instability is RecA-dependent in the (CTG)•(CAG) repeats. Structural studies of the (CCTG)•(CAGG) repeat tract and the flanking sequence, as well as genetic selection assays may reveal the mechanisms responsible for the repeat instability in E. coli, and this may lead to a better understanding of the mechanisms contributing to the human disease state. ^

Application of cell line based genomic predictors to predict response to targeted therapies in breast cancer

Cancer cell lines can be treated with a drug and the molecular comparison of responders and non-responders may yield potential predictors that could be tested in the clinic. It is a bioinformatics challenge to apply the cell line-derived multivariable response predictors to patients who respond to therapy. Using the gene expression data from 23 breast cancer cell lines, I developed three predictors of dasatinib sensitivity by selecting differentially expressed genes and applying different classification algorithms. The performance of these predictors on independent cell lines with known dasatinib response was tested. The predictor based on weighted voting method has the best overall performance. It correctly predicted dasatinib sensitivity in 11 out of 12 (92%) breast and 17 out of 23 (74%) lung cancer cell lines. These predictors were then applied to the gene expression data from 133 breast cancer patients in an attempt to predict how the patients might respond to dasatinib therapy. Two predictors identified 13 patients in common to be dasatinib sensitive. Sixty two percent of these cases are triple negative (ER-negative, HER2-negative and PR-negative) and 76% are double negative. The result is consistent with the findings from other studies, which identified a target population for dasatinib treatment to be triple negative or basal breast cancer subtype. In conclusion, we think that the cell line-derived dasatinib classifiers can be applied to the human patients. ^

FANCM AND FAAP24 MAINTAIN GENOMIC STABILITY THROUGH COOPERATIVE AND UNIQUE FUNCTIONS

Fanconi anemia (FA) is a rare recessive genetic disease with an array of clinical manifestations including multiple congenital abnormalities, progressive bone marrow failure and profound cancer susceptibility. A hallmark of cells derived from FA patients is hypersensitivity to DNA interstrand crosslinking agents such as mitomycin C (MMC) and cisplatin, suggesting that FA- and FA-associated proteins play important roles in protecting cells from DNA interstrand crosslink (ICL) damage. Two genes involved in the FA pathway, FANCM and FAAP24, are of particular interest because they contain DNA interacting domains. However, there are no definitive patient mutations for these two genes, and the resulting lack of human genetic model system renders their functional studies difficult. In this study, I established isogenic human FANCM- and FAAP24-null mutants through homologous replacement-mediated gene targeting in HCT-116 cells, and systematically investigated the functions of FANCM and FAAP24 inchromosome stability, FA pathway activation, DNA damage checkpoint signaling, and ICL repair. I found that the FANCM-/-/FAAP24-/- double mutant was much more sensitive to DNA crosslinking agents than FANCM-/- and FAAP24-/- single mutants, suggesting that FANCM and FAAP24 possess epistatic as well as unique functions in response to ICL damage. I demonstrated that FANCM and FAAP24 coordinately support the activation of FA pathway by promoting chromatin localization of FA core complex and FANCD2 monoubiqutination. They also cooperatively function to suppress sister chromatid exchange and radial chromosome formation, likely by limiting crossovers in recombination repair. In addition, I defined novel non-overlapping functions of FANCM and FAAP24 in response to ICL damage. FAAP24 plays a major role in activating ICL-induced ATR-dependent checkpoint, which is independent of its interaction with FANCM. On the other hand, FANCM promotes recombination-independent ICL repair independently of FAAP24. Mechanistically, FANCM facilitates recruitment of nucleotide excision repair machinery and lesion bypass factors to ICL damage sites through its translocase activity. Collectively, my studies provide mechanistic insights into how genome integrity is both coordinately and independently protected by FANCM and FAAP24.

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

A Genomic Approach to Identify The Notch Pathway as a Putative Tumor Suppressor in Endometrial Cancer

Endometrial cancer is the most common gynecological malignancy and the fourth most frequently diagnosed cancer among women. The molecular changes that distinguish normal endometrium from endometrial carcinoma are not thoroughly understood. Identification of these changes could potentially aid in identifying at-risk women who are especially prone to develop endometrial cancer, such as obese women and women with Lynch Syndrome. A microarray analysis was performed using normal endometrium from thin and obese women and cancerous endometrium from obese women. We validated the differential expression of ten genes whose expression was significantly up-regulated or down-regulated using qRT-PCR. All of the genes had distinct expression levels depending on the endometrial carcinoma histotype. As a result, they could serve as molecular markers to distinguish between normal endometrium and endometrial cancer, as well as between low grade endometrial carcinomas and high grade endometrial carcinomas. Two of the ten genes validated, HEYL and HES1, are down-stream targets of the Notch signaling pathway. HEYL and HES1 were identified by microarray and qRT-PCR to have a significant decrease in expression in endometrial carcinomas compared to normal endometrium. We further analyzed the differential expression of other components of the Notch signaling pathway, Notch4 and Jagged1. They were also identified by qRT-PCR to be significantly down-regulated in endometrial carcinomas compared to normal endometrium. Therefore, we believe the Notch signaling pathway to act as a tumor suppressor in endometrial carcinomas.